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We present a theoretical proposal to couple a single Nitrogen-Vacancy (NV) center to a supercon¬ 
ducting flux qubit (FQ) in the regime where both systems are off resonance. The coupling between 
both quantum devices is achieved through the strong driving of the flux qubit by a classical mi¬ 
crowave field that creates dressed states with an experimentally controlled characteristic frequency. 

We discuss several applications such as controlling the NV center’s state by manipulation of the flux 
qubit, performing the NV center full tomography and using the NV center as a quantum memory. 

The effect of decoherence and its consequences to the proposed applications are also analyzed. Our 
results provide a theoretical framework describing a promising hybrid system for quantum informa¬ 
tion processing, which combines the advantages of fast manipulation and long coherence times. 
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I. INTRODUCTION 

Electronic spins in semiconductors such as NV cen¬ 
ters in diamond or impurities in silicon are characterized 
by their low decoherence rates r s [1-4]. They would be 
excellent candidates for realizing an operating quantum 
processor if only there was a way to reliably entangle dis¬ 
tant spins [5, 6]. Recently, it has been proposed to solve 
this issue by using a superconducting circuit as a quan¬ 
tum bus [7, 8] . This approach would require to reach the 
strong coupling regime where the coupling strength g be¬ 
tween the spin and the circuit is larger than their respec¬ 
tive decoherence rates (g {T s , T c i rcu i t }). To achieve 
this goal, the coupling constant g should be increased 
by several orders of magnitude compared to the current 
state of the art where g is typically of the order of a few 
Hz [9, 10]. 

The coupling constant g can be greatly increased by 
using flux qubits instead of a linear superconducting res¬ 
onator. Indeed, flux qubits [11-18] are characterized by 
a macroscopic permanent current Ip typically of the or¬ 
der of several hundreds of nA. This current flows in the 
loop of the qubit and generates therefore a large magnetic 
dipole which allows an efficient coupling to spins. Bring¬ 
ing the spin at a distance of ~ 20 nm from the flux qubit 
can be achieved by fabricating an ultra-narrow supercon¬ 
ducting wire and aligning it very precisely with the spin 
(see Fig. 1). This would allow increasing the coupling 
by several orders of magnitude and reaching a coupling 
constant g/2ir ~ 100 kHz, a value much larger than re¬ 
cently reported spin decoherence rates [1, 2]. However, 
flux qubits suffer from a severe limitation. Their large 
magnetic dipole makes them very sensitive to flux noise 
which limits their coherence times. To cope with this 
problem, one should tune the flux threading the loop >I> 
precisely at half a flux quantum ($ = $o/2), the so- 


called optimal point of the flux qubit [13]. At this point 
its energy does not depend on the flux at first order and 
is equal to the gap h A. 



FIG. 1. (Color online) Schematic view of a flux qubit coupled 
to a single spin. The flux qubit consists of a superconducting 
loop with a constriction and intersected by four Josephson 
junctions (in red). Three of the junctions are identical while 
the fourth junction is smaller than others. When the flux 
threading the loop <1? is close to half a flux quantum, the 
system behaves as a two level system The Nitrogen Vacancy 
center represented as a golden arrow is situated in the close 
vicinity (r = 15 nm) of the constriction to maximize the cou¬ 
pling between the two systems. The application of an external 
magnetic field R e xt in the direction of the NV axis lifts the 
degeneracy between the ms = ±1 states. 

The coupling of a spin - whose transition energy hco s is 
given by Nature - with a flux qubit is optimal for w s = A. 
Unfortunately, reaching this target is a challenging task. 
The value of A is an exponential function [15] of the pa¬ 
rameters of the junctions which form the qubit and it 
is therefore poorly controlled. It is possible to design a 
tunable flux qubit by replacing one of the junctions by a 
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SQUID [19]. However, coherence times required for sin¬ 
gle spin coupling have not been reported so far in this 
geometry. Applying an external magnetic field to tune 
the spin close to A frequency has also many limitations. 
The critical field of Aluminium is ~ 100 G and there¬ 
fore limits the tunability of the spin frequency should be 
limited to 100-200 MHz. As a consequence the two reso¬ 
nance frequencies w s and A may be highly detuned. Typ¬ 
ically, one can expect that (A — uj s )/2tt ~ 300 — 500 MHz, 
prohibiting the use of the direct resonant interaction to 
entangle the two systems. In this work, we propose a 
theoretical strategy to cope with this issue and couple 
efficiently both systems even with such a detuning. 

II. THE MODEL 
A. Effective Hamiltonian 

We now present our model, starting with NV centers, 
the FQ and the coupling mechanism. NV centers are dia¬ 
mond color centers consisting of a Nitrogen substitutional 
impurity next to a vacancy [20, 21]. In their negatively- 
charged state, they have an electronic ground state with 
a spin S' = 1 which is described by the following Hamil¬ 
tonian (neglecting the effect of strain [22]): 

H s = hDSl + n le S • H ex \ (1) 

where S are the dimensionless spin-1 operators, 
D/(2tt) « 2.88 GHz is the zero-field splitting ground 
state and the rightmost term is the Zeeman interaction 
of the electronic spin of gyromagnetic ratio |7 e |/27r = 28 
GHz/T. The quantization z-axis is set in the Nitrogen 
Vacancy axis. B ext is an external magnetic field applied 
to lift the degeneracy of the {ms = ±1} levels. For 
an antiparallel field, the NV center is well described by a 
two-level system Hamiltonian in the {ms = 0, ms = +1} 
basis: 

H a = hu s a s J 2 , ( 2 ) 

where uj s = D + 7e H ext , & s z = |1)(1| - |0)(0| with |0) and 
|1) corresponding respectively to ms = 0 (lower state) 
and ms = 1. 

The flux qubit Hamiltonian can also be described by a 
two level system [15]: 

H fq = hea f x q /2 + hAa f z q /2 , (3) 

where e = 2/ p ($ — $ 0 / 2 ). The NV is in a position such 
that the field from the flux qubit is normal to the quanti¬ 
zation axis of NV spin S. Therefore the Hamiltonian for 
the two interacting systems can be written as 

H = h,uj s a s z /2 + hea{, q /2 + TiAa{ q /2 + hga x a{ q 7 (4) 

where the fourth term is the Zeeman interaction of the 
electronic spin due to the magnetic field generated by 


the flux qubit. The coupling constant g is given by g ~ 
W I P /(V2- 2nr) where go is the vacuum permeability 
and r the distance between the spin and the constriction 
shown in Fig. 1. By setting r = 15 nm and I p = 500 nA, 
we get g/2ir ~ 100 kHz. In the following, we will assume 
that $ = 4 >q/ 2, in order to benefit from the longest flux- 
qubit coherence times. 

The flux qubit is subjected to an additional microwave 
magnetic field of frequency ui in a direction normal to the 
plane of the persistent currents. The NV spin state has 
a much smaller magnetic dipole and therefore we will 
assume that it is not driven directly by the microwave 
field. Thus the Hamiltonian reads: 

H = Tiu> s a s z /2 + fiAa{ q /2 + hga x a{. q 

+ hHa{ q cos (ut) (5) 

where 12 = IpABo is the Rabi frequency of the flux qubit 
of area A and persistent current Ip driven by the classical 
microwave field B(t) = B 0 cos (cat). 

From (5) we can move to a frame rotating with the 
classical field’s frequency w, which is set to be the same 
as the flux qubit transition frequency - namely uj = A. 
This is accomplished through the unitary operator trans¬ 
formation U s (t) = e — * A *( cr 2+ CT £")/ 2 U r (t), connecting the 
time evolution operators in the Schroedinger and rotat¬ 
ing frame pictures. We define <5 = w s — w = w s — A 
and the raising/lowering operators in the a z basis, a± = 
(a x ± ia y ) /2. By setting 12 ~ |<5|, we identify two terms 
in the resulting Hamiltonian in the rotating frame: 

H r = H 0 + H inU (6) 

with 

H 0 = h S -ai + h^af q ( 7 ) 

Bi n t=ng(a + e + a_e )(o+e +a_e ). 

We can thus consider that the classical drive “dresses” 
the flux qubit [23, 24]. From now on this representation 
will be referred to as the laboratory frame. In order to 
understand the physical meaning of this time-dependent 
interaction Hamiltonian, we plotted in Fig. 2 the time 
evolution of the probability of finding both systems in 
their excited states |1)|1). The figure displays two dis¬ 
tinct time scales: the first one is given by the Rabi oscil¬ 
lations of frequency 12 induced by the classical field. The 
second one comes from the magnetic coupling between 
the flux qubit and the spin. For the numerical simula¬ 
tion we have set oj s = 28800# and A = 25800# so that 
<5 = 3000#, with fl = —5. 

We now move to the interaction picture with respect 
to H 0 in order to get an approximate expression for the 
interaction Hamiltonian. We have: 

Hj = h 9 - |d;e l(2A+5) ‘ [<3-7? _ e if2t («j79 - iaf q ) /2 

+e" int {a{ q + ia{ q ) /2] + a s + e' 5t [&f q - (8) 

e~ int (a{ q + io-79) / 2 + e im (a{ q - ia s y q ) /2] + H.c} . 
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FIG. 2. (Color online) Probability of finding the composite 
system in state 11) 11) as a function of time (in units of 2n/g). 
We numerically solve the Schrodinger equation with Hamilto¬ 
nians Ho and Hint from Eq. 7. We can see the slow dynamics 
induced by the magnetic coupling while the inset displays 
the fast Rabi oscillations occurring with a much shorter time 
scale. The set of parameters is u) s = 28800<?, A = 25800g and 
Q. = -6. 

Let us define = |+)(—| = ( a z — ia y )/2 and <7-, x = 
|-)(+| = (dr + ia y )/2 1 where |±) = (|1) ± |0))/v / 2- We 
recall that we focus on the case with SI ~ (5 in Eq. (8). 
Assuming also the conditions (A,17, <5 g) we can per¬ 
form the so-called rotating wave approximation and are 
led to the effective Hamiltonian: 

H eS (t) = ±7i| (a s + a f ± q x e i(n±s)t + a s _a f T q x e~^ n±s)t ^ , 

(9) 

The effective Hamiltonian (9) describes a closed two 
qubit system that can be diagonalized and studied ana¬ 
lytically. It can lead to a number of applications that will 
be further detailed in the next Sections. We stress that in 
the case of exact matching condition 17 = ±5 the effective 
Hamiltonian becomes time-independent. In the following 

we set fl = — d so that H e g = (a+a+ q x + . 

Adapting the results to SI = 5 is straightforward. 

Supposing that the spin-flux qubit bipartite system is 
initially in a separable pure state, it can be written as: 

|i/>(0)) = (cos0|O) +e* ¥3 sin6 ) |l))(cosa|-t-) + e l<t> sina|—)). 

( 10 ) 

In (10), we dropped subscripts and used a convention 
that will be kept the same throughout this manuscript: 
the spin’s state is expressed in the basis of a s z eigenstates, 
10), 11) and the flux qubit state is expressed in the ba¬ 
sis of a[ q eigenstates, |+), |—). This basis choice is well 
adapted to the coupling induced by (9). Moreover, it 
provides a clear distinction between each party’s states 
without using auxiliary labels. The parameters 6 and ip 
characterizing the spin’s state are unknown and cannot 
be controlled or manipulated without the coupling (9). 


However, the parameters a and (j> can be experimentally 
controlled and the flux qubit can be prepared in an ar¬ 
bitrary initial state by combining classical pulses with 
different phases, intensity and duration [16]. The time 
evolution of the initial state (10) subjected to Hamilto¬ 
nian (9) is: 

| ip(t)} = cos0cosa|O, +) + sin0sina|l, —) 

sin#cosa(—7 sin ^|0, —) + cos yy |1, +)) 

Telcos# sin a (cos yjO, —) — 7 sin yjl,+}). (11) 

In order to test the validity of our model, we compare 
the time evolution induced by the effective Hamiltonian 
(9) to the evolution under the action of the exact Hamil¬ 
tonian in the interaction picture (8) . This can be done by 
comparing the numerical integration of the Schrodinguer 
equation using (8) to the analytical diagonalization of 
(9). We can study as an example the case where 9 = ir/2 
and a = 0. In this case, the initial state is |1)|+) using 
the basis convention. This state is an eigenstate of the 
free Hamiltonian H 0 , so the interaction Hamiltonian is 
solely responsible for the dynamics. The time evolution 
of the population in |1)|+) under either H e g or H lnt is 
shown in Figure 3, with 17 = —S (left) and 17 = —5 + g 
(right). We can see that indeed, both curves display a 
very good overlap. Thus the effective Hamiltonian is a 
good approximation to the coupling between both sys¬ 
tems, even if the resonance is not perfect, i.e. f2 ^ ±<5, 
in this case leading to a reduction of the fringes visibility. 




FIG. 3. (Color online) Probability of finding the composite 
system in state |1)|+) as a function of time (in units of 2n/g). 
In dashed blue we show the exact time evolution of population 
while in continuous red the one obtained using the effective 
Hamiltonian (9). Left: Case of perfect resonance, 17 = —S. 
Right: Non-resonant coupling <5 + 17 = g. We can see that in 
both cases the effective evolution matches the exact one with 
very good accuracy. 

The tuning of the resonance can be done by contin¬ 
uously changing the intensity and phase of the classical 
dressing field until coupling between the NV center and 
the flux qubit is optimal, as in a “quantum radio”. This 
point will be further discussed in the Section III A. 

We should also mention that similar ideas have been 
proposed to couple superconducting qubits with widely 
different frequencies [25, 26]. 
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B. Decoherence 

The dynamics of the coupled spin-flux qubit system 
will be influenced by the presence of noise, that creates 
population losses and dephasing. In order to describe the 
proper dynamical equation for the density matrix of the 
system, we will start from a microscopic description of 
the flux qubit-environment interaction and then derive 
the corresponding Lindblad equation. Given the typical 
decay rates of a flux qubit compared to those of a NV 
center, we will neglect contributions to decoherence in¬ 
duced by the dephasing of the latter. Our microscopic 
model is thus made of two interactions: a dissipative pro¬ 
cess conveyed by the qubit a x operator and a dephasing 
one due to a z . More precisely we have: 

-Hdiss = ^^2'lkH q (Pk + SI) (12) 

k 

and [27] 

H de ph = ^^2 X kV f z q {Ck + 4) ( 13 ) 

k 

where the b^’s an< l 4'^’s are bosonic annihilation (cre¬ 
ation) operators of the modes of the environment. In 
our study we will assume the environment is in thermal 
equilibrium at zero temperature. 

The flux qubit undergoes two independent non- 
dissipative evolutions: the first is due to the microwave 
drive, at frequency Cl, the second to the coupling with 
the NV center at frequency g/2. The latter is the one we 
are interested in while the former is included in the free 
Hamiltonian H 0 . We derive the master equation based 
on these ingredients in the so-called Born-Markov ap¬ 
proximation [28]. The Lindblad equation should model 
the effect of damping and dephasing in the basis of the 
free Hamiltonian eigenstates. Since the free Hamiltonian 
of the flux qubit is proportional to a x , the Lindblad op¬ 
erators will be composed of raising/lowering operators 
in the Pauli-ct x basis. Specifically, we get the following 
dynamical equation for the system density matrix in the 
laboratory frame 

^ = -i[H 0 + H eS ,p] 

+T x C[a f x q ]p + T_C[a f _ q x ]p + T+C[a f £ x ]p (14) 

where C is the Lindblad superoperator C[a\p = a pa ^ — 
^(a^ap+pa^a). Notice that in contrast with the standard 
approach the microwave drive rotates the decoherence 
basis: the ground state is a thermal distribution in the 
{| —), ]-(-)} basis. Even though the environment was as¬ 
sumed to be at zero temperature, the driving held leads 
to an effective heating of the bath. Moreover, the T^’s 
decay rates are related to environmental spectral proper¬ 
ties and thus to experimentally accessible quantities [29]. 
Let us define T\ as the energy relaxation time and T v the 


Rabi oscillations decoherence time which can be mea¬ 
sured through independent experiments. Then we have: 

r * = \t^\ r _ = -t^\ r + = + t~ x ( 15 ) 

One should notice the asymmetry between T + and T_ 
that will have peculiar consequences later on. With the 
dynamical equation (14) we are able to numerically eval¬ 
uate the time evolution of the system, including in our 
simulations realistic values for T\ and T„, experimentally 
determined e.g. in [30]. There we have T) = 10 — 20 
ps and T„ = 10 — 15 ps. In Fig. 4 we show the ef¬ 
fect of decoherence on the reduced flux qubit state ob¬ 
tained after tracing over the spin state. We plot the 
probability of finding the flux qubit in the excited state, 
Pi(t) = Tr spin ((l fq \p(t)\lf q )), where Tr spin denotes the 
partial trace over the spin’s degree of freedom. 



Time (2it/g) 

FIG. 4. (Color online) Probability of finding the flux qubit 
in the excited state as a function of time (in units of 2ir/g), 
for a spin initial state |1). Time evolution is computed in 
the laboratory frame. The plot is at resonance condition, 
decoherence is T\ = 20 ps, T v = 15 ps. Inset: zoom on the 
fast Rabi oscillations for very short times. 

We will now show how this formalism can be used in 
fundamental quantum information protocols. 

III. APPLICATIONS IN QUANTUM 
INFORMATION PROTOCOLS 

We will see in the following how a proper choice of 
the parameters a and <fi can lead to a number of ap¬ 
plications relying on the coupling between the two-level 
systems, such as manipulation and measurement of the 
spin’s state. 

A. Detecting spin-qubit coupling 

The first step in view of exploring the flux qubit-single 
spin coupling is to detect the coupling itself. Here we 
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study a protocol enabling this in presence of decoherence, 
which introduces noise and loss of information about the 
flux qubit’s state. 

We first provide the exact calculation neglecting dis¬ 
sipation and restricting to zero detuning, i.e. = —6, 
focusing on the interaction picture. We define P{{t) = 
Tr spin ((l/q|/3*(i)|l/q)), where p l (t) is the density matrix 
of the whole system computed in the frame rotating with 
Hq. Based on Eq. 11 we have: 


P[{t) 


1 

2 


^1 + sin 2a cos cj> cos 



(16) 


We recall that the main goal of the present protocol is to 
detect the existence of the coupling. In this regard, the 
key point of the dynamics of the population P[{t) is that 
it does not depend on the NV center state, which is a 
priori unknown, but rather on the flux qubit state which 
is controllable. This feature will not be preserved when 
including decoherence to the dynamics. Hence we can set 
a = 7t/ 4 (i.e. starting with the flux qubit in the excited 
state |1)) to get the highest visibility. In that case P[(t) 
reads: 


pm 


1 

2 


gt 

1 + cos 


(17) 


This equation is valid provided the flux qubit and the 
spin are in perfect resonance Q = —S. Furthermore, the 
resonance condition ensures that P\{t) corresponds to 
the enveloppe of the fast oscillations when moving back 
to the laboratory frame (see Fig 5, top left). Eq. (17) 
corresponds to an oscillation with unit visibility. The 
fast oscillations at frequency fl come from the classical 
driving field which induces the free Hamiltonian of Eq. 7. 
They are therefore of no interest regarding the coupling 
with the spin. 

In practice, exact resonance must be determined ex¬ 
perimentally. In this respect, driving the flux qubit off- 
resonance reduces the visibility of the coupling oscilla¬ 
tions, as can be seen in Fig. 5. So one practical way to 
find out the resonance is by modifying II and searching 
for high visibility oscillations. It is possible to gradually 
vary H until the slowly oscillating terms reach a popula¬ 
tion inversion dynamics - practically cross the P\ = 1/2 
horizontal line, which uniquely characterizes the reso¬ 
nance condition. Fig. 5 shows the time evolution of Pi{t), 
defined in section IIB, for different values of the detun¬ 
ing. It is clearly possible to identify when the resonance 
condition is matched. Moreover, the amplitude of the os¬ 
cillations characterizing the coupling are very sensitive to 
the detuning - typically for values as low as g/2 - which 
guarantees the precision of this protocol. 

However, decoherence may significantly damp the os¬ 
cillations and destroy the expected signal related to the 
population measurement of state |1). Furthermore, as 
will be shown later, the dynamics gets to depend on the 
spin’s initial state. In order to characterize the effects 
of decoherence we studied the behavior of P\(t) in the 
presence of decoherence for different values of the system 
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FIG. 5. (Color online) Probability of finding the reduced flux 
qubit density matrix in the excited state |1) as a function of 
time (in units of 27r/gi). Decoherence here is not taken into 
account. From top left and increasing clockwise, values of the 
detuning are: + 8 = 0, g/2, g and 2 g. 


coherence times and detuning (see Fig. 6). We set the 
initial state of the spin to the ground state |0) and based 
the numerical analysis on the Lindblad equation derived 
in Eq. 14. It is clear from Fig. 6 that even when includ¬ 
ing decoherence, the coupling of the flux qubit with the 
spin can be identified through the measurement of os¬ 
cillations in the excitation probability of the flux qubit. 
The sensitivity to the detuning is sustained, which en¬ 
ables searching for the resonance condition by gradually 
varying Q. 

One may ask how low the coherence times can be in or¬ 
der to see the slow oscillations due to the coupling with 
the spin. The coupling should be strong enough with 
respect to decoherence to allow for the experimental ob¬ 
servation of at least one oscillation. In Fig. 7, we show 
the oscillations at resonance condition for decreasing val¬ 
ues of Ti, while T„ is set to 15 ps. The initial state is 
|0 S ) 11 / q )- Surprisingly enough, the figures reveal that the 
coupling can be identified for rather low values of T) such 
as 2 ps, that is gT\ ~ 0.2. 

An interesting feature of the coupling also appears 
when dissipation is considered: depending on the spin 
initial state, it can either protect the flux qubit from the 
effects of decoherence or enhance them. This effect is 
illustrated in Fig. 8 where two different initial states of 
the spin are considered, leading to different damping be¬ 
haviors. This is a consequence of the specific form of the 
master equation Eq. (14), which exhibits an asymmetry 
between the dissipation terms corresponding to <r+ :X and 
1 7 - tX . We could indeed verify that the dependence on 
the spin’s initial state vanishes when artificially impos¬ 
ing T+ = T_. 
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FIG. 6. (Color online) Probability of finding the reduced flux 
qubit density matrix in the excited state |1) as a function 
of time (in units of 2n/g). Decoherence is modeled through 
Ti = 20 fis, T v = 15 fis. Black crosses show the enveloppe of 
the fast oscillations in the absence of coupling g = 0. From 
top left, and increasing clockwise, values of the detuning are 
: SI + 5 = 0, g/ 2, g and 2 g. 
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FIG. 7. (Color online) Probability of finding the reduced flux 
qubit density matrix in the excited state |1) as a function of 
time (in units of 2iv/g). T„ = 15 gs while Ti is 20, 10, 5 
and 2 gs (from top left, decreasing clockwise). Black crosses 
show the enveloppe of the fast oscillations in the absence of 
coupling <? = 0. 

B. Spin state initialization 

A first application of the coupling described in (9), that 
is also the first of the DiVincenzo [31] criteria to define a 
qubit, is the spin state initialization. By setting a = 0 in 
(10), we have that the initial state of the hybrid system 
is given by 

(18) 



Time (2jt/g) Time (2n/g) 

FIG. 8. (Color online) Probability of finding the reduced flux 
qubit density matrix in the excited state |1) as a function 
of time (in units of 2n/g). Decoherence is modeled through 
Ti = 20 gs, Tv = 15 gs. Black crosses show the enveloppe of 
the fast oscillations in the absence of coupling g = 0. Left: 
spin initial state is |0). Right: spin initial state is |1). 


From the time evolution induced by the coupling between 
the spin and the flux qubit, after t = w/g , this state will 
be set to state |0). From this point on, the NV center 
can either be manipulated and set to an arbitrary state, 
or can be used as a quantum memory, as developed in 
the next subsection. 

We now analyze the impact of decoherence on the ini¬ 
tialization protocol. The results are shown in the left 
column of Table I. Even though we consider the NV cen¬ 
ter itself to be perfectly isolated from the environment, 
the purity is reduced by the interaction with the flux 
qubit and leads to the preparation of non-pure state in¬ 
stead of the target one |0). In Table I we call fidelity 
Trfi ux ((0 s |/5(7r/g)|0 s )), where Trfl ux denotes the partial 
trace over the flux qubit’s degree of freedom. In the 
presence of high decoherence rates, one way to improve 
the state initialization is by repeating the protocol. Af¬ 
ter a spin-flux qubit interaction time of t = 7 r/g, one 
measures the flux qubit state disregarding the results 
of the measurement. Mathematically, this measurement 
projects the flux qubit and destroys entanglement with 
the spin, which corresponds to tracing over the measure¬ 
ment eigenbasis. Then one may again prepare the flux 
qubit in |+). We then let the spin and the flux qubit 
interact for the same duration. The repetition of the pro¬ 
tocol improves the state initialization procedure, until it 
reaches an asymptotic value independent of the initial 
spin state - right column of Table I. 

The values of the fidelity between the final reduced 
density matrix of the spin and the projector on its ground 
state - see Table I - were numerically calculated based 
on the dynamical evolution equation (14). In order to 
compute those values, we repeated the protocol described 
in the former paragraph and then averaged the fidelity 
we obtained over 500 Haar-random initial spin states. 


C. Quantum memory and spin state manipulation 

We have shown in Section III B how to initially prepare 
the spin in state |0). From this initialized state, we show 


|^(0)) = (cos0|O) + e lip sin0|l))|+). 
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Typical times 
(Ti,T„) (in ns) 

Average fidelity 
of the protocol 

Average fidelity 
after 5 iterations 

(10,10) 

0.92 

0.95 

(20,15) 

0.96 

0.97 


TABLE I. Average fidelities for the initialization protocol. 
Details on how they were computed can be found in the main 
text. 

now how the spin can be used as a quantum memory, 
encoding in a long lifetime quantum system the state of 
the flux qubit. We will also see how to adapt this strategy 
to manipulate the spin’s state, realizing arbitrary single¬ 
qubit rotations. 

We start with the quantum memory protocol. If the 
initial spin state is |0), we have that 9 = 0 in (10) and the 
flux qubit is prepared in an arbitrary state that should 
be perfectly transferred to the spin. 

IV’(O)) = |0)(cosa|+) + e*^ sina|—)). (19) 

Using (11), we see that after a time t = ir/g with the 
initial condition (19), the SC quantum state will be com¬ 
pletely transferred to the spin state. Thus, the latter can 
play the role of a quantum memory, since it has a longer 
coherence time than the flux qubit. Figure 9 shows the 
fidelity of the final spin state after this protocol with re¬ 
spect to the initial flux qubit state when decoherence is 
also considered. The computed fidelity does not depend 
on the phase <f>, because the decoherence process itself 
is (/(-independent. The maximum fidelity is obtained for 
a = 0, that is for an initial state |0) |+) which is an eigen¬ 
state of the coupling Hamiltonian of eigenvalue 0. The 
minimum is for a = 7t/ 2 which implies a complete ex¬ 
citation transfer and is therefore the most likely to be 
impacted by decoherence. 

The strategy above enables the spin state manipulation 
as well. Suppose that one wants to apply a given rotation 
to an arbitrary spin state, that can be expressed as in 
(18). Any rotation can be associated with the following 
transformation of the basis states of the NV center: 

|0) = 1 6 ) —> cos /?|0) + e iX sin /3| 1) 

|1) = |1) —> cos/?|l) — e _iX sin/3|0). (20) 

This transformation can be achieved by the following pro¬ 
tocol: in a first step, the NV center state is reinitialized, 
as described in IIIB. Thus, the state of the final coupled 
system composed by the spin and flux qubit is given by: 

It k) = |O)(cos0|+) + e lv sin0|—)), (21) 

be., the state of the spin is completely transferred to the 
flux qubit. After state transfer, the coupling between the 
spin and the flux qubit can be stopped by turning off the 
intense dressing classical pulse. The flux qubit can then 
be manipulated by another classical field as follows: 

|+) —> cos/3|+) + e* x sin/3| —) 

|—) —> cos/3| —) — e“* x sin/3|+). (22) 



FIG. 9. (Color online) Fidelity of the NV center reduced 
density matrix after the protocol with the initial flux qubit 
state (a) as a function of a and </> (b) as a function of a 
only, (a) shows that the fidelity of the protocol is actually 
(/(-independent. Values of decoherence for the two plots in (b) 
are the same as in Table I. 


Finally, the dressing microwave field can be turned on, 
coupling again the spin to the flux qubit. After t = 7r /g, 
we obtain the final state: 

\ip r ) = cos0|O) + e llp sin0|l), (23) 

which corresponds exactly to the realization of an arbi¬ 
trary rotation to the initial state (18). 


D. NV center state tomography 

The previously discussed strategies can also be used 
for realizing the full spin state tomography. Since we 
have shown that it is possible to transfer the - unknown 
- state of the NV-center to the flux qubit, one can, after 
this operation, simply realize the full tomography of the 
flux qubit. 

Suppose now the initial state is such that a = 0 and 
<j> = 0, so the flux qubit is in |+). Then for a pulse dura- 

















tion t = 7r /g, the state becomes, according to Eq. (11): 

| i/j(t)) = |O)(cos0|+) + ie llp sin0|—)) (24) 

The unknown state of the NV-center has been transferred 
to the accessible flux qubit. Full tomography of the latter 
yields perfect knowledge about the initial state of the 
NV-center. 

We have also studied the fidelity of this protocol in 
presence of decoherence, finding that it displays the 
same behavior as the one of the initialization protocol. 


IV. CONCLUSIONS 

We have shown that a hybrid system, composed of 
a directly inaccessible spin and a superconducting flux 
qubit can be effectively coupled by driving the flux qubit 
with an intense classical microwave field even in the limit 
where both two-level systems have far off resonant char¬ 


acteristic frequencies. The coupling is created by the 
dressing of the flux qubit by the classical field. Such 
dressing leads to the possibility of tuning the dressed 
eigenvalues to the frequency difference between the two 
quantum devices, a process that can be described by an 
effective Hamiltonian. The possibility of coupling such 
devices, that present complementary advantages with re¬ 
spect to quantum information processing, leads to a num¬ 
ber of applications, discussed in the present paper. We 
have developed protocols to manipulate the spin state, 
use it as a quantum memory and realize its full tomog¬ 
raphy. In all protocols, a detailed study of the effects 
of decoherence in the dressed system was included, es¬ 
tablishing limits on the expected fidelities according to 
decoherence rates compatible to the state of the art. Our 
results serve as a roadmap to promising experiments us¬ 
ing hybrid quantum devices. An interesting perspective 
is studying how the flux qubit can serve as a data bus, 
intermediating the coupling between the spin and the 
quantum field of a resonator, both in the strong and the 
ultra-strong limits [32, 33]. 
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